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ABSTRACT 

We present the results on millimeter interferometric observations of four lu- 
minous infrared galaxies (LIRGs), Arp 220, Mrk 231, IRAS 08572+3915, and 
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VV 114, and one Wolf-Rayet galaxy, He 2-10, using the Nobeyama Millime- 
ter Array (NMA). Both the HCN(l-O) and HCO+(1-0) molecular lines were 
observed simultaneously and their brightness-temperature ratios were derived. 
High-quality infrared L-band (2.8-4.1 /im) spectra were also obtained for the 
four LIRGs to better constrain their energy sources deeply buried in dust and 
molecular gas. When combined with other LIRGs we have previously observed 
with NMA, the final sample comprised nine LIRGs (12 LIRGs' nuclei) with avail- 
able interferometric HCN(l-O) and HCO+(1-0) data — sufficient to investigate 
the overall trend in comparison with known AGNs and starburst galaxies. We 
found that LIRGs with luminous buried AGN signatures at other wavelengths 
tend to show high HCN(1-0)/HCO^(1-0) brightness-temperature ratios as seen 
in AGN-dominated galaxies, while the Wolf-Rayet galaxy He 2-10 displays a 
small ratio. An enhanced HCN abundance in the interstellar gas surrounding a 
strongly X-ray-emitting AGN, as predicted by some chemical calculations, is a 
natural explanation of our results. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: ISM — radio 
lines: galaxies — galaxies: individual (Arp 220, Mrk 231, IRAS 08572-1-3915, VV 
114 and He 2-10) 



1. Introduction 



It has recently been found that supermassive blackholes (SMBH) with masses of Msmbh 
> lO^M© are ubi quitously present at t he centers of spheroidal components (bulges and 
elliptical galaxies) (IMagorrian et al. I ll998h . When active mass accretion onto a SMBH occurs, 
it is observed as an active galactic nucleus (AGN). The curr ent AGN picture postulates that 
dust is present in torus geometry around the central SMBH (jAntonuccilll993l ) . The direction 
perpendicular to the torus is mostly transparent to the bulk of the AGN's ionizing UV - 
soft X-ray photons. Gas clouds at 10-1000 pc distances from the AGN along the torus axis, 
photo-ionized by the central AGN's radiation, produce the so-called narrow-line regions 
(NLRs; Robson 1996). Since the spectral shapes of the ionizing photons differ between 
AGNs and stars, emission-line ratios from NLRs differ from gas clouds photo-ionized in star- 
forming galaxies. Thus, AGNs surrounded by torus-shaped d ust with well developed NLR s 
are distinguishable from star-forming galaxies relatively easily (jVeilleux fc Osterbrocklll987l ). 



Because dust at the central <10 pc of galaxies has an angular momentum, the spatial 
distribution of the dust is likely to be axisymmetrical; dust column density is higher in one 
direction (the torus direction) and lower in another direction (the torus axis). When the 
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nuclear concentration of dust becomes large, even the torus axis direction can be opaque 
to the bulk of the AGN's UV - soft X-ray radiation, and the radiation can be blocked at 
the inner part (<10 pc), producing virtually no NLRs. Such buried AGNs are no longer 
detectable as long as emission lines originating in the NLRs are sought. Because a simple 
interpretation of the cosmic X-ray backgroun d spectrum suggests that most AGNs in the 
universe are buried (IFabian fc Iwasawal 119991 ) . it is important to establish a technique for 
finding buried AGNs. In a buried AGN, almost all of the energetic UV - soft X-ray radi- 
ation is absorbed by the surrounding dust and is re-emitted as infrared dust emission. For 
this reason, luminous buried AGNs are expected to be strong infrared emitters. Luminous 
infrared galaxies (LIRGs), which radiate strong infrared emission with Li^ > IO^^Lq and 



contain highly concentrated nuclear dust (ISanders fc Mirabellll996l ). are probable objects in 
which luminous buried AGNs reside. 

Since star-forming activity can also produce infrared emission, it is necessary to dis- 
tinguish the origin of the infrared luminosities of LIRGs if we are to study the luminous 
buried AGN population. An effective method is to search for the presence of strong X-ray 
emission because an AGN intrinsica lly produces much stronger X-ray emission compared to 
star-forming activity (IRobson 



U nfortun ately, buried AGNs generally contain a large 



amount of nuclear dust (llmanishi et al.ll2006a[). and most of them may be Compton-thick 



(Nh > 10 cm ) sources (iMaiolino et al.l l2003l ) . X-ray observations at E > 10 keV are 



needed to directly detect Compton-t hick X-ray emission from luminous buried AGNs with a 
small scattered/reflected component (IFabian et al.ll2002l ): however, these are applicable only 
to a few very nearby bright LIRGs due to a current sensitivity limit. 

Although a systematic buried AGN search through direct X-ray observations is dif- 
ficult, we can investigate the presence of a luminous buried AGN through the chemical 
effects of the intrinsically strong X-ray emission on the surrounding interstellar medium. 
Around a strongly X-ray-emitting buried AGN surrounded by dense gas and dust, the so- 
called X-ray dissociation regions (XDR; Maloney et al. 1996) should develop, in contrast 
to photo-dissociation regions (PDRs) usually seen in strongly UV-emitting star-forming re- 
gions. Detection of XDR signatures can thus be a use ful tool for locatin g a luminous burie d 
AGN. Among the several proposed probes of XDRs (lUsero et al.ll2004l : lAalto et al.l 120071 ). 
we foc us on the observationaUy derived method based on HCN(l-O) and HCO"'"(1-0) lin e 
ratios (IKohno et al.l l200ll : iKohnd l2005l : llmanishi et al.l l2006bl : iGracia-Carpio et al.l |2006| ) , 
because t heoretical predictiori of line fluxes from PDRs and XDRs requires many free pa- 
rameters (lYamada et al.l 120071 ) and could still be uncertain, given that the r nost important 
param eter — detailed properties of the clumpy structure of molecular gas (ISolomon et al. 
19871 ) — is totally unknown from observations. 
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Here, we present the results of our millimeter interferometric HCN(l-O) and HC0"''(1- 
0) observations using the Nobeyama Millimeter Array (NMA) of nearby LIRGs whose en- 
ergy sources have previously been investigated at other wavelengths. Compared to pre- 
yious HCN(l-O) and HCO+fl-O) observations of LIRGs using single-dish radio telescopes 
( iGao fc Solomon! l2004l : iGracia-Carpio et al.l |2006| ) , our interferometric maps have the im- 
portant advantage that spatially resolved HCN(l-O) and HCO"'"(1 -0) data are obtainable i n 
LIRGs, which often show disturbed, multiple nuclei morphologies (ISanders &: Mirabellll996l ). 
New ground-based infrared L-band (2.8-4.1 /zm) spectra of higher quality than older data 
were also taken for individual main nuclei of these millimeter-observed LIRGs to better con- 
strain the nature of their obscured energy sources. We adopted Hq = 75 km s~^ Mpc~^, Qm 
= 0.3, and = 0.7 throughout this paper. 



2. Targets 



The target LIRGs are selected on the basis of their proximity {z < 0.06) to cover the 
redshifted HCN(l-O) and IICO"^(1-0) emission lines with the NMA receiving systems, and 
their expected high fluxes of HCN(l-O) and HCO+(1-0). The four LIRGs, Arp 220, IRAS 
08572-f3915, Mrk 231, and VV 114, were observed. Table 1 summarizes the infrared emission 
properties of these LIRGs. An angular scale of 1" corresponds to a physical size of 0.3-1 kpc 
at redshift = 0.018-0.058. 

Arp 220 (z = 0.018) is one of the best-studied nearby ultraluminous infrared galaxies 

(ULIRGs; LiR > lO^^L©; Sanders et al. 1988a ). It consists of two n uclei, Arp 220 E 

and W, with a separation of ~1" (jScoville et al.ll2000l : ISoifer et al.ll2000l ). and the optical 



spectrum of the combined emission is clas sified as a LINER ( i.e., non-Seyfert; no obvious 
AGN signatures in the optical spectrum) (IVeilleux et al.lll999l ). The presenc e of starburst 
activ i ty (= active star-formati on) is evident in various observational data (iGenzel et al. 



19981 : llmanishi &: Dudleyl l2000l ). but the detected starburst is energetically insufficient to 



fully account for the observed infrared lumino sity of Arp 220 quantitatively (llmanishi et al. 



Arp 220's nuclear core (ISpoon et al 



2004 



2006al : lArmus et al.ll2007l: llmanishi et al.ll2007f). requiring the energy sou rce deeply buried in 



Gonzalez- Alfonso et al.l 120041 ). The presence of a 



luminous buried AGN at the core 



radio ( Shioya et al 



2001 



las long bee n unclear even through detailed observ ations 



in the X-rav (jlwasawa et al.ll2005l). infrar ed (IGenzel et al.l 1 19981 : lArmus et al.l 120071 ) , and 



Parra et al.ll2007l ) frequencies. However, the recent discovery of a 



compact energy source with a very high emission surface brightness has provided the first 
strong signatures of a luminous buried AGN in Arp 220W (jPownes fc Eckartll2007l ). 



Mrk 231 (z = 0.042) is a well studied, single nucleus ULIRG and is the most luminous 
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object (LiR ~ lO^^' ^LfT^) in the local un iverse. It is an optically known AGN classified as 
a Seyfer t 1 galaxy (IVeilleux et al.lll999[). but it displays a bsorption features in the infrared 
spectra ( Roche et al. 19831 : Armus et al. 2007 ) and X-rays JSraito et al "2004h possibly orig- 
inating in the broad absorption line (BAL) clouds (IBoksenberg et al.iil977i ). The presence 
of a very luminous and probably energetically dominant AGN is revealed by various obser- 
vations (e.g., Soifer et al. 2000; Imanishi & Dudley 2000; Braito et al. 2004; Imanishi et al. 
2006a). 



IRAS 08572+3915 (z = 0.058) i s a ULIRG (L 
and SE) with a separation of ~5" (IScoville et al. 



2000 



10 ^^'^L(7)) corisistin g of two nuclei (NW 
Kim et al.l[2002[). The NW nucleus 



(IRAS 08572+3915NW) is believed to be energetically dominant fISoifer et al.ll2000l) and is 
therefore our primary interest. It is classified optically as a LINER (jVeilleux et al.lll999l ). 
It is one of the strong buried AGN candidates because (1) the energy source is suggested 
to be very compac t and is more centrally concentrated than the nuclear dust, as expected 
for a buried AGN JPudlev fc Wvnn-WiUiam]ll997l : ISoifer et al.lboool : llmanishi et aPboOGal . 
20071 ). and (2) the polyc yclic aromatic hydrocarbons (PAH) emission (the starburst indica- 



Armus et al. 2007; Imanishi et al 



tor) is extreraely weak ( Imanishi &: Dudlevi l2000l ; llmanishi et al.l l2006al ; ISpoon et al.l 12006 



200 



3). It is thus a particularly interesting target for in 



vestigating whether the HCN(l-O) to HCO"'"(1-0) brightness-temperatu re ratio is sim ilar to 
those found in AGN-dominated nuclei or starburst-dominated galaxies (jKohndl2005l ). 



V V 114 (z = 0.020 ) is a double-nuclei LIRG (Lir ~ with a separation of 



15" (iKnop et al.lll994j ). Both nuclei are classified optically as HII regions (IVeilleux et al. 



19951 ). In the short optical wavelength, the western nucleus (VV 114W) is brighter than 
the eastern nucleus (VV 114E), but VV 114E b ecomes more important with increasin g 
wavelength and is prominent in the infrared band (IKnop et al.lll994l : iLe Floc'h et al.ll2002l ). 
probably dominating the infrared luminosity of the VV 114 merging system. The VV 114E 
nucleus has a secondary peak ~1'.'5 southwest of the VV 114E peak (VV 114Esw; Soifer et 
al. 2001). The presence of a luminous buried AGN i s suggested in VV 114 E from the strong 
featureless mid-infrared 15 fim continuum emission (ILe Floc'h et al.ll2002l ). 



In Kohno (2005), HCN(l-O) to IICO^(l-O) brightness-temperature ratios are studied in 
AGNs and typical starburst gal axies such as M82. In the nuclei of LIRGs, where molecular 
gas is very highly concentrated ( Sanders fc Mirabe3ll996l ). putative star formation could be 
extreme in that young, massive, hot stars are more predominant than a normal starburst. To 
determine if such an extreme starburst shows a different HCN(1-0)/HCO^(1-0) brightness- 
temperature ratio from a normal starburst, a nearby well stud ied Wolf-Rayet g alaxy. He 
2-10 (z = 0.003), dominated by massive, hot Wolf-Rayet stars (I Allen et al .11 19761 ) was also 
studied (Table 1). 
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3. Observations and Data Reduction 

3.1. Millimeter interferometry 

Interferometric observations of HCN(l-O) (Arest = 3.3848 mm and z/rest = 88.632 GHz) 
and HCO+(1-0) (Ajcst = 3.3637 mm or i/^cst = 89.188 GHz) lines were undertaken with the 
Nobeyama Millimeter Array (NMA) and the RAINBOW interferometer at the Nobeyama 
Radio Observatory (NRO) between 2005 January and 2007 April. Table 2 summarizes the 
detailed observation log. The NMA consists of six 10-m antennas and observations were 
undertaken using the AB (the longest baseline was 351 m), C (163 m), and D (82 m) con- 
figurations. The RAINBOW interferometer is a seven-element combined array that includes 
the NRO 45-m telescope in addition to the six 10-m antennas (NMA) . The RAINBOW ob- 
servations were scheduled when the NMA array was set at the AB configuration and the 
longest baseline was 410 m. In the 3-mm wavelength range, the sensitivity of the RAIN- 
BOW interferometer is about twice better than that of the NMA array only because the total 
aperture size increases by a factor of four with the inclusion of the NRO 45-m telescope. 



The backend was the Ultra- Wide-Band Correlator (UWBC) flOkumura et al.ll2000l ) con- 
figured to cover 1024 MHz with 128 channels at 8-MHz resolution. The central frequency 
for each source (Table 2) was set to cover both the redshifted HCN(l-O) and HCO+(1-0) 
lines simultaneously. A bandwidth of 1024 MHz corresponds to ~3500 km s~^ at ~ 84-89 
GHz. The field of view at this frequency is ~26" (full-width at half-maximum; FWHM) and 
~77" (FWHM) for the RAINBOW and the NMA array, respectively. The Banning window 
function was applied to reduce side lobes in the spectra; thus, the effective resolution was 
widened to 16 MHz or 55 km s~^ at ~ 84-89 GHz. Since the declinations of VV 114 and 
He 2-10 are low (< —15°), an observable time period (elevation > 20°) per each day from 
the NMA site is short, requiring many observing days. 



The UVPROC-II package developed at NRO fiTsutsumi et al.lll997f ) and the AIPS pack- 
age of the National Radio Astronomy Observatory were used for standard data reduction. 
Corrections for the antenna baselines, band-pass properties, and the time variation in the 
visibility amplitude and phase were applied to all of the data (Table 2). A fraction of the 
data with large phase scatter due to poor millimeter seeing was removed from our anal- 
ysis. After clipping a small fraction of data of unusually high amplitude, the data were 
Fourier-transformed using a natural UV weighting. The fiux calibration was made using 
observations of Uranus or appropriate quasars (Table 2) whose fiux levels relative to Uranus 
or Neptune had been measured at least every month in the NMA observing seasons. A 
conventional CLEAN method was applied to deconvolve the synthesized beam pattern. The 
primary beam pattern of the NMA antenna was corrected for LIRGs with very extended 
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spatial structures (VV 114 and He 2-10). Table 3 summarizes the total net on-source in- 
tegration times and synthesized beam patterns. Absolute positional uncertainties of the 
NMA/RAINBOW maps were estimated to be much less than 1". 



3.2. Infrared 2.8—4.1 /im (L-band) spectroscopy 



We conducted ground-based infrared 2.8-4.1 /im (L-band) slit spectroscopy on the main 
nuclei of the four LIRGs (Arp 220, Mrk 231, IRAS 08572+3915, and VV 114) to better 
understand the nature of their obscured energy sources based on the equivalent width of 
the 3.3 /im PAH emission and the optical depths of absorption features at Arcst ~ 3.05 /im 
(in the rest frame) by ice-covered dust grains and at Arsst ~ 3. 4 /im by bare carbonaceous 
dust grains (jimanishi fc Dudleyll200d : llmanishi fc Maloneyll2003l ). A normal starburst galaxy 
should always sh ow large equivalent-width PAH emission, while a p ure AGN produces a PAH- 
free continuum (jimanishi fc Dudley! l2000l : llmanishi et al.l l2006al ). For absorption features, 
optical depths have upper limits in a normal starburst in which stellar energy sources and 
dust are spatially well mixed, while they can be arbitrarily lar ge in a buried AGN with a 



more centrally concen trated energy source geometry than dust (jimanishi fc Maloneyl 12003 



Imanishi et al.ll2006al ). For these reasons, these PAH emission and dust-absorption features 
can be used to distinguish the dust-obscured energy sources of LIRGs. 

Ground-based infrared spectra taken with large (>4-m) telescopes are superior in spatial 
resolution to those obtained by space satellites with small diameters (<1 m), thus enabling 
us to obtain spatially resolved spectra of double nuclei with a small separation (<2 arc- 
sec). For example, Arp 220 and VV 114E have double nuclei with separations of 1 to 
1.5 arcsec, which can be resolved only with the ground-based spectra. For Mrk 231 and 
IRAS 0857 2+3915, although infrared L-band spectra taken with <4m telescopes have been 
published ( jimanishi et al.l Il998l : llmanishi fc Dudleyl l2000l : llmanishi et al.l l2006a] ) , a higher- 
quality L-band spectrum obtained with a larger, 8-m-class telescope may be able to provide 
new constraints on their energy sources. We thus obtained infrared L-band spect ra of Arp 



220, Mrk 231, IRAS 08572+3915, and VV 114E using the Subaru 8.2-m telescope (live et al. 
20041 ) atop Mauna Kea, Hawaii. 



We used the IRCS infrared spectrograph (jKobayashi et al.ll2000l ) at the Nasmyth focus 
of the Subaru telescope to obtain the new infrared L-band spectra of these LIRGs. Table 4 
tabulates the observation log. The sky was clear during the observations. The seeing in the 
K-hemd measured in images taken before the L-band spectroscopy was 0'.'4-0'.'8 in FWHM. 
A O'.'G- or O'.'9-wide slit and the L-grism were used with a 52-mas-pixel scale. The achievable 
spectral resolution is R ~ 140-200 at A ~ 3.5 /im. When the seeing size was sufficiently 
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small, we basically used the 0'.'6 slit to search for AGN signatures at a n uclear core based on 
the PAH equivalent width and the strengths of dust abs orption features (llmanishi fc Dudley 
2OOOI : llmanishi &: Maloneyl l2003l : llmanishi et al.ll2006al ) with a reduced contamination from 
extended starburst components. However, for IRAS 08572+3915, the 0'.'9 slit was employed 
because the seeing in the i^-band worsened during the observation. Thus, PAH emission from 
an extended (>0'.'6 or >0'.'9) starburst component is not covered; thus, PAH flux or luminosity 
is not considered in our discussion. For Arp 220, Mrk 231, and IRAS 08572+3915, nuclear 
PAH fluxes (or luminosities) within the centra l 1-2 arcsec can better be investigated from 
previously obtained spectra using a wider slit (llmanishi fc Dudleyll2000l ). The precipitable 
water was low, 1-2 mm for the 2006 July observing run, and <1 mm for the 2007 January 
run. We employed a standard telescope nodding technique (ABBA pattern) with a throw of 
5 to 7 arcsec along the slit to subtract background emission. We used the optical guider of 
the Subaru telescope to monitor the telescope tracking. Exposure time was 1.0-1.5 sec, and 
30-60 coadds were made at each nod position. 

For both observing runs, F- or G-type main sequence stars (Table 4) were observed 
as standard stars, with a mean airmass difference of <0.1 to the individual LIRG nuclei 
to correct for the transmission of Earth's atmosphere and to provide flux calibration. The 
L-band magnitudes of the standard stars were estimated from their l/-band (A = 0.6 fim) 
magni tudes, adopting the V — L colors appropriate to the stellar types of individual standard 
stars jTokunagalboOoi ). 

Standard data analysis procedures were employed using IRAF0. Initially, frames taken 
with an A (or B) beam were subtracted from frames subsequently taken with a B (or A) 
beam, and the resulting subtracted frames were added and divided by a spectroscopic flat 
image. Bad pixels and pixels hit by cosmic rays were then replaced with the interpolated 
values of the surrounding pixels. Finally, the spectra of the LIRGs' nuclei and standard stars 
were extracted by integrating signals over 0'.'8-l'.'8, depending on the actual signal profiles. To 
create a combined spectrum of Arp 220E+W, we summed signals over 2'.'5-2'.'9 along the slit. 
Wavelength calibration was performed based on the wavelength-dependent transmission of 
Earth's atmosphere. The spectra of the LIRGs' nuclei were divided by the observed spectra 
of standard stars and multiplied by the spectra of blackbodies with temperatures appropriate 
to individual standard stars (Table 4). 



Flux calibration was performed based on signals of the LIRGs and standard stars de- 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc. (AURA), under cooperative agreement with the 
National Science Foundation. 
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tected inside our slit spectra. Seeing sizes in the K-band (and L-band) were always smaller 
than the employed slit widths, and good telescope tracking performance of Subaru was es- 
tablished. We thus expected possible slit loss for the compact (<0'.'6) emission to be minimal, 
say <30%. In fact, we divided the whole data set into subgroups and compared their flux 
levels, which agreed within 20% for all LIRGs. This level of possible flux uncertainty (<30%) 
for the compact emission will not seriously affect our main conclusions. To reduce scatter of 
the infrared spectra, appropriate binning of spectral elements was performed, particularly at 
Aobs < 3.3 /im in the observed frame, where the scatter of data points is larger due to poorer 
Earth atmospheric transmission than at Aobs > 3.4 /zm. 



4. Results 

4.1. Millimeter interferometric data 

For Arp 220 and Mrk 231, millimeter spectra at HCN(l-O) or HCO"'"(1-0) emission peak 
positions show that the flux levels between these lines are substantially above the zero level, 
indicating that continuum emission is clearly present. We combined data points that are 
unaffected by these lines and made interferometric maps of the continuum emission. Figure 
1 presents the contours of the continuum emission. The continuum emission is clearly (>9cr) 
detected in Arp 220 and Mrk 231. We barely see the continuum emission signature for VV 
114, IRAS 08572+3915, and He 2-10, but its detection significance is <3ct. Table 5 presents 
the estimated continuum levels. 

Figure 2 shows the integrated intensity maps of the HCN(l-O) and HCO^(l-O) emission 
of the five observed sources (Arp 220, Mrk 231, IRAS 08572+3915, VV 114, and He 2-10). 
The continuum emission is subtracted for Arp 220, Mrk 231, and VV 114, but not for IRAS 
08572+3915 and He 2-10 because of possible large ambiguities. 

Since molecular gas is highly concentrated at the nuclei of LIRGs, t he importance of 
high- density molecular gas (uh > 10^ cm~^) is expected to increase there fiGao &: Solomon 



20041 ). Such molecular gas is better probed with HCN(l-O) and HCO^(l-O) rather than the 
widely used CO (1-0) because the dipole moments of HCN and HCO+ (/i > 3 debye) are 
much larger than CO (/i ~ 0.1 debye; Botschwina et al. 1993; Millar et al. 1997). Hence, 
the maps in Figure 2 reflect the spatial distribution of this important dense molecular gas. 
A spatially unresolved emission peak is clearly seen for Mrk 231 at the nuclear position. For 
IRAS 08572+3915, the spatially unresolved HCN(l-O) emission peaks at the energetically 
dominant NW nucleus, as expected (see § 2, paragraph 4). For He 2-10, H CO"'"(1-0) emission 



displaying no significant spatial extent is found close to the CO (1-0) peak (IKobulnicky et al. 
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19951 ). A spatially extended emission component is detected for Arp 220 and VV 114. Figures 
3 and 4 display the channel maps around the HCN(l-O) and HCO^(l-O) lines of Arp 220 
and VV 114, respectively. 

For Arp 220, the HCN(l-O) and HCO+(1-0) emission peaks reside in between the 
double nuclei in the integrated intensity maps of Figure 2. Nevertheless, the channel map 
in Figure 3 clearly shows that Arp 220E becomes more prominent with increasing velocity 
(upper-left, red component) and Arp 220W does with decreasing velocity (lower-right, blue 
component). Our HCN(l-O) and HCO+(1-0) channel maps thus confirm that the Arp 220E 



nucleus is farther redshif t ed than Arp 220W, as was previously noted (jLarkin et al.l 11995 
Downes &: Solomon 1998 : Taniguchi &: Shioya 1998 : Sakamoto et al. 19991 ). 



In the case of VV 114, the HCN(l-O) the emission peak is close to VV 114E, while 
that of HCO"''(1-0) is strong between the VV 114E and W nuclei (Figure 2), similar to 
CO (1-0), (2-1), and (3-2) emissi 



(Yun et al. 


1994: 


lono et al. 


2004) 



butions of HCN(l-O) and HCO+(1-0) are quite different, possibly because (a) HCN(l-O) 
and HCO+(1-0) probe molecular gas with different density (iGreve et al.ll2006l ). due to their 
slightly different critical densities □, or (b) HCO ^ (1-0) may selectivel y trace shock regions. 



due to enhanced HCO^ abundance by shocks (jPickinson et al.lll980l ). Further details will 
be discussed in a future paper (Tamura et al., in preparation). No significant HCN(l-O) or 
HCO"^(1-0) emission is recognizable in VV 114W. 

Figure 5 shows spectra around the HCN(l-O) and HCO^(l-O) lines for Arp 220, Mrk 
231, IRAS 08572+3915, and VV 114. For Mrk 231 and IRAS 08572+3915, spectra are 
extracted at the nuclear peak position. For Arp 220, the spectra are extracted at both the 
Arp 220E and W positions, and double-peaked HCN(l-O) and HCO^(l-O) emission is seen 
at both positions. Although the double-peaked profile suggests that emission from Arp 220E 
and W are not clearly spatially resolved in a map, the redshifted (lower frequency) component 
is stronger at the Arp 220E position than at Arp 220W, as expected from the channel map 
(Figure 3). We can thus separate the emission from Arp 220E and W spectroscopically, with 
the red and blue components coming from Arp 220E and W, respectively. For VV 114, both 
the HCN(l-O) and HCO'^(l-O) emission have spatial structures, and their peak positions 



The critic al density (iic) is, b y definition, proportional to Einstein's A coefficient (A) divided by the 



collision rates ( Schoier et al.ll2005r) . The A-coefficient A is oc /i^ x , where [i is the dipole moment and v 
is the frequency of the molecular line. The dipole moments (/i) of HCN (/i ~ 3 debye; Millar et al. 1997) 
and HCO+ (/x 4 debye; Botschwina et al. 1993) are similar to each other. The frequencies of HCN(l-O) 
(88.632 GHz) and HCO+(1-0) (89.188 GHz) are almost identical. Thus, the A -coef5cient is s i milar between 
HCN(l-O) and HCO+(1-0). However, using the collision rates estimated by ISchoier et all (j2005r ). the Uc 
value of HCN(l-O) is a factor of -5 larger than that of HCO+(1-0). 
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are not the same. We thus extract spectra at four different positions: (E-1) the HCN(l-O) 
peak at the eastern side of VV 114E in Figure 2; (E-2) the HCN(l-O) peak at the western 
side of VV 114E in Figure 2; (3) the eastern HCO+(1-0) peak between VV 114E and W in 
Figure 2; and (4) the western HCO"'"(1-0) peak between VV 114E and W in Figure 2. 

Figure 6 presents Gaussian fits to the detected HCN(l-O) and HCO"'"(1-0) hues. The 
central velocity and line width of the Gaussian fit are determined independently for the 
HCN(l-O) and HCO+(1-0) lines. For IRAS 08572+3915, a double Gaussian fit was at- 
tempted because the lines seem to be double-peaked. Table 6 summarizes the Gaussian 
fitting results. For Mrk 231 and IRAS 08572+3915 (the LIRGs dominated by the single pri- 
mary nucleus), the HCN(l-O) peak velocity is simi l ar to that of CO(l -O) in previous reports 
JPownes fc SolomonlEggsl : ISrvant fc Scovillelll99d : lEvans et al. l l2002h . For Arp 220E and W 
nuclei, the peak velocities of HCN(l- O) and HCO+(1-0) agree with that of CO (2-1) within 
50 km s~^ (IDownes fc Solomonlll998l ). 



For Arp 220, Mrk 231, IRAS 08572+3915, and VV 114, the integrated fluxes of HCN(1- 
0) and HCO'^(l-O) at each peak position estimated from the Gaussian fits in the spectra 
are summarized in Table 7. For He 2-10, the fluxes are derived from the peak contours 
of the integrated intensity maps (Figure 2) because the spectrum is noisy. The HCN(1- 
0)/HCO"'"(1-0) brightness-temperature ratios (oc x flux density) are also shown in Table 
7. In Mrk 231 and IRAS 08572+3915, both the HCN(l-O) and HCO+(1-0) emission show 
a spatially unresolved single peak at the main nuclear position. The integrated HCN(l-O) 
and HCO^(l-O) fluxes are also estimated from the peak contours of the integrated intensity 
maps in Figure 2. The estimated values are consistent with those based on the Gaussian 
fits for both HCN(l-O) and HCO+(1-0) emission within 15% in Mrk 231 and 30% in IRAS 
08572+3915 if the continuum in Figure 6 is assumed. 

Arp 220, VV 114, and He 2-10 display spatially extended HCN(l-O) and HCO+(1-0) 
emission. Their total fluxes are shown in Table 8. For Arp 220 and Mrk 231, HCN(l-O) 
and HCO~'"(1-0) fluxes have previously been measured using single-dish radio telescopes. 
Table 9 compares their fluxes with our measurements, which are similar, suggesting that our 
interferometric data recover the bulk of HCN(l-O) and HCO^(l-O) emission. 



4.2. Infrared spectra 



4-2.1. PAH emission 



Figure 7 presents the flux-calibrated infrared 2.8-4.1 /xm (L-band) spectra of Arp 220, 
Mrk 231, IRAS 08572+3915NW, and VV 114E. For Arp 220 and VV 114E, our ground-based 
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spectra provide spatially resolved spectra of double nuclei with small (1-2 arcsec) separation, 
which is not possible for space-based infrared spectra (see §3.2). Thus, our L-band spectra 
make it possible to discuss the energy sources of individual nuclei sepa rately for Arp 220 E 
and W JScoville et al.ll200ol : ISoifer et al.l[200o[ ). and VV 114E and Egw dSoifer et al.ll200lh . 



The 3.3 /im PAH emission is clearly seen in the spectra of the observed LIRGs, except 
for IRAS 08572+3915NW, indicating that a detectable amount of starburst activity is surely 
present. Table 10 summarizes the strength of the 3.3 fim PAH emission feature. The ob- 
served 3.3 fj,m PAH luminosities (L3 3pah) roughly trace the absolute magnitudes of modestly 
obscured (Ay < 15 mag) starburst activities covered inside our slit spectra (the Lir of such 
a starburst is ~10^ x Ls.spah; Mouri et al. 1990; Imanishi 2002). 

For Arp 220E+W, Mrk 231, and IRAS 08572+3915NW, the equivalent widths of the 3.3 
fim PA H emission (EW.s .spah) are similar to previous es timates based on old, lower-quality 
spectra (jimanishi fc Dudley! l200d : llmanishi et al.ll2006al ). The EWs.spah values of Mrk 231, 
IRAS 08572-I-3915NW, and W 114E are more than a factor of ~5 smaller than the typical 
value found in starburst-dominated galaxies (EWs spah ~ 100 nm; Imanishi & Dudley 2000), 
indicating that PAH-free continua from the putative luminous AGNs contribute irnportantly 
to th eir observed L-band fluxes ( llmanishi fc Dudleyl l2000l : llmanishi et al.l l2006al : llmanishi 
20061 ). For Arp 220E, W, and VV 114Esw5 the EW3.3PAH values are as large as seen in 
normal starburst galaxies, so that no explicit signs of PAH-free continua from buried AGNs 
are detected in this wavelength range. 

For Mrk 231 and IRAS 08572-1-3915, the ULIRGs dominated by compact nuclear emis- 
sion, the L-band continuum flux levels in our slit s pectra are similar , within 0.1 mag, to 
the photometric measurements using a 2'.'5 aperture (IZhou et al.lll993l ). For Arp 220E-I-W, 
the L-band continuum flux level is a factor of ~2 smaller than previously obtained spectra 



using a 1^.^2 apertu re (llmanishi fc Dudleyl I2OOOI ) and photometric data using a 2'.'5 aperture 
( IZhou et al.lll993l ). The small flux is possibly due to the small slit width (0'.'6) in the new 
spectra and/or possible L-band flux ambiguity in the adopted standard star, HR 5728 (Table 
4). 



4-2.2. Absorption features 



IRAS 08572-1-39 15NW (Figure 7) shows a strong absorption feature at Aobs = 3.6 /im 
in the observed frame or Arest = 3.4 /um in the rest frame due to bare carbonaceous dust 
grains. Its optical depth is estimated to be Txa 0.8, which is similar to previous estimates 
( llmanishi &: Dudleyll2000l : llmanishi et al.ll2006al ). The optical depth is larger than the upper 
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limit produced with the natural geometry of a normal starburst (stellar energy sources and 
dust spatially well mixed; r3.4 < 0.2), but requires buried AG N-type geometry in which 



the energy source is more centrally concentrated than the dust (jimanishi fc Maloneyl 12003 



Imanishi et al.ll2006al ). supporting the buried- AGN classification as previously indicated from 



the low EW3.3PAH value (<3 nm; Table 10). 

In the high-quality IRAS 08572+3915NW spectrum in Figure 7, the broad absorption 
feature centered at Aobs = 3.2 /im or Arest = 3.0-3.05 /xm is seen for the first time. In 
a buried AGN and starburst composite galaxy, the buried AGN emission is more highly 
obscured, because the central compact buried AGN should be geometrically surrounded 
by the starburst. Hence, the contribution from buried AGN emission to an observed fiux 
becomes higher with increasing wavelength, from the iC-band (2-2.5 ^m) to L-band (2.8- 
4.1 /Um), due to smaller dust extinction. Additionally, blackbody radiation from stars (T 
> 4000K) generally shows a steeply decreasing fiux from the K- to L-band. For these 
reasons, it could happen that an observed i^-band continuum emission largely comes from 
foreground stellar emission, in addition to buried AGN emission, while an L-band continuum, 
particularly at a longer wavelength part, is dominated by buried AGN emission. If this is 
the case for IRAS 08572 -I-3915NW and if the stellar emission's tail extends from the K-band 
to the shorter part of the L-band (Aobs < 3.0 /im), then an absorption-like feature at Arest ~ 
3.05 /xm could be reproduced. However, the K-band spectrum of IRAS 08572-1-3915 comes 
mostly from buried-AGN-heated hot dust emission, because of the e xtremely red /T-band 



continuum and no st ellar CO absorption features at Arest ~ 2.3 /im (iGoldader et al.lll995 



Imanishi et al.ll2006al ). Thus, this scenario seems implausible. We ascribe the absorption- like 
feature to the 3.05 /xm absorption by ice-covered dust grains and estimate its optical depth 
to be Ts.i ~ 0.3. 

The spectrum of VV 114E also displays the 3.05 /xm ice-absorption feature with an 
optical depth of 1 ~ 0.5. Its optical depth is again larger than the threshold explained 
by the mixed dust/source geometry of a normal starburst (T3.1 < 0.3; Imanishi & Maloney 
2003), suggesting the presence of a luminous buried AGN, as inferred from the low EW3.3PAH 
value (~20 nm; Table 10). 

For the infrared L-band spectra of the remaining LIRG nuclei (Arp 220E, W, Mrk 231, 
and VV 114Esw), no clear absorption features are recognizable. The absence of absorption 
features in Mrk 231 (the low EW3.3PAH object; Table 10) implies that the AGN is only weakly 
obscured, while those for Arp 220E, W, and VV 114Egw are explained by the predominant 
contribution from starburst emission to the ob served L-band fiuxes. For Arp 220W, the 



nucleus with luminous buried AGN signatures (IDownes fc EckartI 120071 ) . we see no explicit 



AGN sign in the infrared L-band spectrum. The emission from the putative buried AGN 
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may be so highly flux-attenuated that its contribution to an observed L-band flux must be 
insignificant, compared to foreground weakly-obscured starburst emission. 



Discussion 



5.1. Comparison of HCN(1— 0)/HCO+(1— 0) brightness-temperature ratios 

with other galaxies 

Figure 8 shows the HCN(1-0)/HCO+(1-0) and HCN(1-0)/CO(1-0) brightness-temperature 
ratios for Arp 220, Mrk 231, IRAS 08572+3915, VV 114 and He 2-10. Previously ob tained 



data points of nearby LIRGs { 



manishi et al. 2004, 2006b: Imanishi fc Nakanishi 2006). star- 



bursts, and Seyfert galaxies (IKohnd l2005l ) are also plotted. As stated by llmanishi et al. 



(l2006bl ). we mainly use the HCN(1-0)/HCO^(1-0) ratios (ordinate) in our discussions for 
the following reasons. First, since the HCN(l-O) and HCO+(1-0) lines are observed simul- 
taneously with the same array configuration of NMA/RAINBOW, their beam patterns are 
virtually identical. We can thus be confident that the same regions are probed with both 
lines. Second, both HCN(l-O) and HCO'''(1-0) fluxes are measured at the same time with 
the same receiver and same correlator unit so that the possible absolute flux calibration un- 
certainty of interferometric data does not propagate to the uncertainty in the ratio, which is 
dominated by statistical noise and fitting errors (see Figure 6). Therefore, the derivation of 
the HCN(1-0)/HCO~^(1-0) brightness-temperature ratio is straightforward. Neither of the 
first and second arguments hold for the HCN(1-0)/CO(1-0) brightness-temperature ratios in 
the abscissa of Figure 8, which complicates their interpretation. LIRGs with luminous AGN 
signatures at other wavelengths (Arp 220W, Mrk 231, IRAS 08572+3915) tend to show high 
HCN(1-0)/HCO"^(1-0) brightness-temperature ratios as seen in AGN-dominated galaxies. 
The ratio is also high in VV 114E-1, which is close to the putative AGN location at VV 
114E, despite a slight positional deviation. The Wolf-Rayet galaxy He 2-10 shows a very 
low HCN(1-0)/HCO^(1-0) brightness-temperature ratio. Some individual examples of the 
nearby Galactic high-mass-star-fo rming, Hll-region cores show high HCN (1-0)/HCO"'"(1- 
0) brightness-temperature ratios (ITurner &: Thaddeud 119771 : iPirogovl Il999l ). However, for 
external galaxies, the ensemble of HII regions, molecular gas, and PDRs will be observed 
as a whole. The low HCN(1-0)/HCO+(1-0) brightness-temperature ratio in He 2-10 sug- 
gests that even an extreme starburst dominated by young hot stars cannot easily produce a 
high HCN(1-0)/HCO^(1-0) brightness-temperature ratio when observed at a large physical 
scale. The high ratios observed in LIRGs are likely to have another cause than the extreme 
starburst scenario. 
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5.2. Interpretations of HCN(1— 0)/HCO+(1— 0) brightness-temperature ratios 



A natural explanation for the high HCN(l-O) /HCO"*"(1-0) brightness-temperature ratios 
in buried AGN candidates is the increasing abundance of HCN due t o strong X-ray radiation 



from the AGN , as p redicted by some simple chemical calculations (iLepp fc Dalgarndll996 
Lintott fc VitilboOGi ). When both the HCN(l-O) and HCO+(1-0) emission are optically 
thin, the HCN(1-0)/HCO^(1-0) brightness- temperature ratio will increase linearly with 



increasing HCN abundance relative to HCO"'". Ilmanishi et al.l (j2006bl ) and Figure 9 show that 
even if the emission is moderately optically thick, the increasing HCN abundance results in 
a larger HCN(1-0)/HCO"'"(1-0) brightness-temperature ratio (see also Knudsen et al. 2007). 
Since the HCN(1-0)/HCO+(1-0) brightness-temperature ratios in Figure 8 derived from our 
NMA/RAINBOW data are only toward the cores where putative AGNs are expected to be 
present, with a minimum contamination from extended starburst emission outside the beam 
size, the potential AGN's X-ray chemistry, if any, may better be reflected than single-dish 
telescope's measurements. 

The critical density (uc) HCN(l-O) is a factor of ~5 larger than that of HC0^(1- 
0) (§ 4.1). When the fraction of high density gas, above the HCN(l-O) critical density, 
increases, the HCN(1-0)/HCO"'"(1-0) brightness-temperature ratio could be high compared 
to low density gas. Sinc e LIRGs have highly concentrated molecular gas in their nuclei 



(jSanders fc Mirabel! 1996 ) and the density may be larger than the normal starburst galaxies 
studied by Kohnol (2005), the large HCN(1-0)/HCO+(1-0) brightness-temperature ratios 
could originate simply in the increased molecular gas density. However, in this scenario, 
the HCN(1-0)/CO(1-0) brightness-temperature ratios should be correspondingly high (see 
§4.1). Although HCN(1-0)/CO(1-0) brightness-temperature ratios in the abscissa of Figure 
8 exhibit some ambiguities (§ 5.1), no such expected trend is observed in Figure 8 (LIRGs 
are not distributed in the right side), which does not support this high gas-density scenario. 
The expected high gas density in the nuclei of LIRGs and non-high HCN(1-0)/CO(1-0) 
brightness-temperature ratios can be reconciled if the higher average gas density in certain 
areas results from a larger volume- filling factor of the clump (Figure 9), rather than the 
higher density of each clump. 



Garcia-Burillo et al.l (120061 ) found a very large HCN(5-4) /HCO^(5-4) brightness-temperature 



ratio in a luminous AGN with strong mid-infrared 12 /im continuum emission (APM 08279-1-5255), 
and suggested that infrared radiative pumping, rather than coUisional excitation, can produce 
a high HCN/HCO"'' brightness-temperature ratio. In general, an AGN is a more luminous 
mid-infrared 12 /im continuum emitter than a starburst, because 12 /im hot dust emission 
is stronger in the former. Hence, the infrared pumping may work more effectively in LIRGs 
with luminous buried AGNs than those without, which could also explain the high HCN(1- 
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0)/HCO^(1-0) brightness-temperature ratios seen in luminous buried AGN candidates. 



Papadopoulod (120071 ) proposed that the HCN(1-0)/HCO"'"(1-0) brightness-temperature 



ratio could increase with increasing turbulence of molecular gas. In this model, it is expected 
that the ratio is higher with larger line widths of HCN(l-O) and HCO^(l-O). Figure 10 
compares the line width with the ratio, but we see no clear positive trend. 



Gracia-Carpio et all (120061 ) measured the HCN(l-O) and HCO+(1-0) fluxes of LIRGs 



using a single-dish telescope with >25-arcsec beam size, and found a positive correlation 
between the HCN(1-0)/HCO+(1-0) brightness-temperature ratio and infrared luminosity. 
Figure 11 shows the HCN(1-0)/HCO+(1-0) brightness-temperature ratio only for the core 
emission inside 1.5-10 arcsec beam sizes as a function of infrared luminosity from the whole 
galactic region. Among LIRGs plotted in Figure 11, the two LIRGs d ominated by a single 



main nucleus, Mrk 231 and Mrk 273, were plotted in the figure of iGracia-Carpio et al. 



(120061). The two LIRGs with multiple main nuclei, Arp 220 and Arp 299, were also studied 
by iGracia-Carpio et al.l (120061 ) . but our interferometric data resolved the nuclei, providing 
essentially new plots. Other LIRGs in Figure 11 are new. A positive correlation is also 
recognizable here with a different sample. If our AGN's X-ray interpretation for the origin 
of high HCN(1-0)/HCO^(1-0) brightness-temperature ratios is correct, then an increasing 
importance of AGNs at the cores of LIRGs with higher infrared luminosities is suggested. 
Investigating the spatial variation of the ratios using a better spatial- resolution map (i.e., 
separating a true compact core and surrounding extended regions) may help to clarify the 
origin of the ratios. Detailed theoretical calculations that solve chemical reaction networks 
in a self-consistent way and realistically account for the clumpy structure of molecular gas 
are also desirable for interpreting our observational results. 



6. Summary 

We presented the results of millimeter interferometric IICN(l-O) and HCO^(l-O) ob- 
servations of nearby LIRGs using the NM A/RAINBOW array. When combined with other 
LIRGs we have previously observed, ours is the largest interferometric HCN(l-O) and HC0"'"(1- 
0) survey of LIRGs. Most of the observed LIRGs are selected based on the presence of lu- 
minous AGN signatures at other wavelengths. From the interferometric data, we extracted 
the HCN(l-O) and HCO+(1-0) flux at the core positions, where the putative AGNs are 
expected to be present. For observationally investigating the possible chemical effects of 
the putative strongly X-ray-emitting AGN to the surrounding interstellar gas, our data set 
is much more advantageous than single-dish telescope's data because the contamination 
from extended starburst emission can be reduced. We derived the HCN(1-0)/HCO"^(1-0) 
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brightness-temperature ratios of LIRGs and compared them to the ratios found in AGNs 
and starburst galaxies. 

We reached the following main conclusions. 

1. LIRGs with detected luminous AGN signatures through observations at other wave- 
lengths generally show high HCN(1-0)/HCO''"(1-0) brightness-temperature ratios and 
are distributed in the region occupied by AGNs. 

2. When combined with HCN(1-0)/CO(1-0) brightness-temperature ratios and the hne 
widths of HCN(l-O) and HCO+(1-0), the high HCN(1-0)/HCO+(1-0) brightness- 
temperature ratios observed in these LIRGs are naturally explained by an HCN abun- 
dance enhancement as suggested by some simple theoretical models. 

3. We found a positive correlation between the HCN(l-O) /HCO+(1-0) brightness-temperature 
ratios only toward the LIRGs' cores and infrared luminosities of LIRGs, suggesting that 
AGN activity is more important in galaxies with higher infrared luminosities. 
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Table 1. Detailed information about the observed LIRGs and the Wolf-Rayet galaxy He 

2-10 



Object 


Redshift 


/12 


/25 


/so 


/lOO 


log LiR (log Lir/Lq) 


Far-infrared 






(Jy) 


(Jy) 


(Jy) 


(Jy) 


(ergs s"-') 


Color 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Arp 220 


0.018 


0.48 


7.92 


103.33 


112.40 


45.7 (12.1) 


0.08 (cool) 


Mrk 231 


0.042 


1.87 


8.66 


31.99 


30.29 


46.1 (12.5) 


0.27 (warm) 


IRAS 08572+3915 


0.058 


0.32 


1.70 


7.43 


4.59 


45.7 (12.1) 


0.23 (warm) 


VV 114 (IC 1623) 


0.020 


0.68 


3.57 


22.58 


30.37 


45.2 (11.6) 


0.16 (cool) 


He 2-10 


0.003 


1.09 


6.51 


24.08 


26.40 


43.5 (9.9) 


0.27 (warm) 



Note. — Column (1): Object name. Column (2): Redshift. Columns (3)-(6): fi2, f25, feOi and fioo are IRAS 
FSC (or PSO) fluxes at 12, 25, 60, and 100 fim, respectively. Column (7): Decimal logarithm of the infrared 
(8-1000 Atm) luminosity in ergs s'^ calc ulated as follows: Lir = 2.1 X 10^^ X D(Mpc)2 x (13.48 X /12 + 
5.16 X /25 + 2.58 X feo + /loo) ergs s"! llSanders fc Mirabel 199^. The values in parentheses are the decimal 
logarithms of the infrared luminosities in units of solar luminosities. Colu mn (8): IRAS 25 ^m to 60 fim flux 
ratio (f25/f6o)- LIRGs with f25/f60 < (>) 0.2 are classified as cool (warm) jSanders et al.iri988bl ). 

''Distance of 10.5 Mpc ||Tu11vII1988|) is adopted. 



Table 2. Observing log 



Object 
(1) 


Array 
configuration 
(2) 


Observing Date 
(UT) 
(3) 


Central 
frequency 
(4) 


Phase 
calibrator 
(5) 


Band-pass 
calibrator 

(6) 


Flux 
calibrator 
(7) 


Arp 220 


RAINBOW 
AB 


2006 Jan 26, 28, 30 

2007 Jan 22, 23, 24, 25 


87.34 


1538+149 


3C 273 


1741-038 


Mrk 231 


RAINBOW 


2006 Jan 25, 26, 28 


85.31 


1418+546 


3C 273 


1741-038 


IRAS 08572+3915 


RAINBOW 
AB 
C 
D 


2005 Jan 28 

2006 Feb 12 

2005 Mar 13, 14, 25, Dec 11, 13 

2006 Mar 4, 5 


84.04 


0923+392 


3C 84, 3C 273 


3C 84 


VV 114 


AB 
C 
D 


2007 Jan 29, 30, Feb 6, 7, 8 
2007 Mar 29, 30 

2006 Nov 28, 29, 30, Dec 1 

2007 Apr 26, 27 


87.16 


0048-097 


3C 84, 3C 454.3 


Uranus 


He 2-10 


AB 
C 
D 


2007 Feb 6 

2007 Mar 28, 29, 30 

2006 Dec 13, 2007 Apr 26, 27 


88.64 


0834-201 


3C 84, 3C 279 


3C 84 



Note. — Column (1): Object name. Column (2): NMA array configuration. Column (3): Observing date in UT. Column (4): 
Central frequency used for the observations. Column (5); Object name used as a phase calibrator. Column (6): Object name used as 
a band-pass calibrator. Column (7): Object name used as a fiux calibrator. 
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Table 3. Parameters of final NMA/RAINBOW map 



Object 


On source 


Beam size 


Position ang 


;Ie of 




integration (hr) 


(arcsec x arcsec) 


the beam 


n 


(1) 


(2) 


(3) 


(4) 




Arp 220 


24 


1.8 X 1.5 


-55.1 




Mrk 231 


4 


2.1 X 1.8 


-34.8 




IRAS 08572+3915 


33 


4.7 X 3.5 


-50.7 




VV 114 


18 


7.5 X 5.5 


-13.2 




He 2-10 


13 


10.8 X 5.5 


-10.1 





Note. — Col.(l): Object name. Col.(2): Net on-source integration time in 
hours. Col.(3): Beam size in arcsec x arcsec. Col.(4): Position angle of the beam 
pattern. It is 0° for the north-south direction, and increases counterclockwise 
on the sky plane. 



Table 4. Subaru IRCS observing log 



Object 


Date 


Integration 


Slit width 


P.A. 




Standard Stars 






(UT) 


(Min) 


(") 


n 


Name 


L-mag 


Type 


Toff (K) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Arp 220 


2006 July 17 


24 


0.6 


87 


HR 5728 


4.5 


G3V 


5800 




2007 January 14 


20 


0.6 


87 


HR 5728 


4.5 


G3V 


5800 


Mrk 231 


2007 January 14 


16 


0.6 


90 


HR 4767 


4.8 


F8-G0V 


6000 


IRAS 08572+3915NW 


2007 January 15 


24 


0.9 


90 


HR 3451 


5.0 


F7V 


6240 


VV 114E 


2006 July 19 


24 


0.6 


60 


HR 210 


4.0 


G3V 


5800 



Note. — Column (1): Object name. Col. (2): Observing date in UT. Col. (3): Net on-source integration time in min. Col. 
(4): Slit width in arcsec. Col. (5): Position angle of the slit. 0° corresponds to the north-south direction. Position angle increases 
counterclockwise on the sky plane. Col. (6): Standard star name. Col. (7): Adopted L-band magnitude. Col. (8): Stellaj: 
spectral type. Col. (9): Effective temperature. 
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Table 5. Continuum emission 



Object 


Flux 


Frequency 




(mJy) 


(GHz) 


Arp 220 


30 


-87 


Mrk 231 


35 


-85 


IRAS 08572+3915NW 


<3 (< 3a) 


~84 


VV 114 


<4 (< 3a) 


-87 


He 2-10 


<6 (< 3a) 


-88.5 



Table 6. Gaussian fitting parameters of HCN(l-O) and HCO"*"(1-0) emission lines 



Object 




LSR velocity 




FWHM 






(km s-l) 




(km 


s-l) 




HON (1-0) 


HCO+(1-0) 


CO(l-O) or 00(2-1) 


HON(l-O) 


HCO+(1-0) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Arp 220E 


5640 


5680 


5650 


240 


190 


Arp 220W 


5310 


5290 


5340 


320 


270 


Mrk 231 


12660 


12690 


12650-12660 


260 


290 


IRAS 08572+3915NW 


17110 ^ + 17540 


17250 ^ -1- 17900 


17490 


420 =1 + 300 


320 -f 180 


VV 114E-1 


6090 






170 




VV 114E-2 


6070 


6060 




190 


260 


VV 114-3 


6060 


6030 




290 


360 


VV 114-4 


5910 


5950 




250 


310 



Note. — Col.(l): Object name. Col.(2): LSR velocity {vopt = (^ - 1) X c} of the HCN(l-O) emission peak. Col.(3): 
LSR velocity o f the HCO^(l— 0) emission peak. OoI.(4): LSR velocity of the 00(1—0 ) or 00(2— 1) emission peak, taken from 
the hterature l|Brvant fc Scoville|[l99i : lOownes fc Solomonlll998l : lEvans et aljbooj ). Col.(5): Line width of the HON (1-0) 
emission at FWHM. Ool.(6): Line width of the HCO+(1-0) emission at FWHM. 

''A double Gaussian fit. 
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Table 7. Flux of HCN(l-O), HCO+(1-0), and CO(l-O) emission at each peak position 



Nucleus 


HCN(l-O) 


HCO+(1-0) 


HCN(1-0)/HCO+(1-0) 


CO(l-O) 


Reference 




(Jy km s~^) 


(Jy km s~^) 




(Jy km s~^) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Arp 220E 


11.9 


5.3 


2.3 


30 ^ 


A 


Arp 220W 


17.8 


6.6 


2.7 


47 ^ 


A 


Mrk 231 


10.7 


6.3 


1.7 


68 


B 


IRAS 08572+3915NW 


2.3 (1.4 + 0.9) ^ 


1.3 (0.7 + 0.6) ^ 


1.8 


11.8 


C 


VV 114E-1 


2.7 


<1.7 


>1.6 


<70 


D 


VV 114E-2 


1.9 


4.0 


0.5 






VV 114-3 


3.2 


5.5 


0.6 






VV 114-4 


1.3 


6.0 


0.2 






He 2-10 


6.1 


11 d 


<0.6 


30 


E 



Note. — Col.(l): Object name. Col. (2): Integrated HCN(l-O) intensity at an emission peak. Col. (3): Integrated 
HCO+(l— 0) intensity at an emission peak. Col. (4): HCN(1-0)/HCO+ (1-0) ratio in brightness temperature (oc X flux 
density) at an emission peak. The ratio is not affected by possible absolute flux uncertainties in the NMA/RAINBOW 
data (see §5.1). Col. (5): Integrated CO(l-O) intensity a t an emission pe a k take n from the literature. Col. (6): Ref erences 
for interferometric CO( l-O) flux measur ement s. fA): ISakamoto et aL lll999h . (B); iDownes fc SolomonI lll998h . (C): 



lEvans et all ||2002|) (D): llono et"al] ||2004| ). (E): iKobulnickv et all \l993t ) 



^^00(2-1) emission is detected at individual cores, b ut CO(l-O) emission is not jPownes fc Solomonll998|). We estimat e 
the C0(1— 0) flux at the core from the CO(2-l) flux jSakamoto et al.l ll999"). using the scaling llRiechers et alj|2006h . 
based on the assumption that emission is thermalized and that both CO(l-O) and CO(2-l) have the same brightness 
temperatures. 

''Blue (former) and red (latter) components from the double Gaussian fits (see §4.1). 
^^00(1-0) fiux is estimated from CO(3-2) flux, by assuming the scaling. 

"^Flux is derived from the peak contour of the integrated intensity map. Possible continuum emission is not subtracted, 
because it is difficult to estimate in a reliable way due to the large scatter of the spectrum. 

''The ratio is an upper limit because possible continuum emission is not subtracted and HCN(l-O) emission is weaker 
than HCO+(l— 0) emission. 



Table 8. Total flux of HCN(l-O) and HCO"'"(1-0) emission for sources that show 

spatially-extended components 



Nucleus 


HCN(l-O) 


HCO+(1-0) 




(Jy km s^-"-) 


(Jy km s"-"-) 


(1) 


(2) 


(3) 


Arp 220 


34 


18 


VV 114 


7 


19 


He 2-10 


6.5 


11 



Note. — Col.(l): Object name. Col.(2): 
Total HCN(l-O) flux. Col.(3): Total 
HCO+(1-0) flux. 
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Table 9. Comparison of HCN(l-O) and HCO^(l-O) fluxes between our measurements and 

tfiose in the literature 



Nucleus 


Line 


Flux 


Reference 




(1) 


(2) 


(3) 


(4) 




Arp 220 


HCN(l-O) 


34=^ 


This work 








45 


Evans et al. (2006) 








37 


Solomon et al. ( 1992' 


1 






35=^ 


Radford et al. (1991") 






HCO+(1-0) 


IS'' 


This work 








23 


Gracia-Caroio et al. I 


(2006) 


Mrk 231 


HCN(l-O) 


ll^' 


This work 








14 


Solomon et al. ( 1992' 


1 




HCO+(1-0) 


6=^ 


This work 








8 


Gracia-Caroio et al. ( 


f 20061 


Note. 


— Col.(l): Object name. Col.(2): HCN(l-O) or 


HCO+(1-0) line. Col.(3): 


Flux in 


[Jy km s-l]. Col.(4): 


Refer- 



ence. 



''Measurements are made based on interferometric data. 



Table 10. Strength of the 3.3 /xm PAH emission feature 



Object 


fs.SPAH 


L3.3PAH 


L3.3PAH/L1R 


rest EW3.3PAH 




(xlO^-'-* ergs cm~^) 


(xlO^'^ergs s-i) 


(xlO-3) 


(nm) 


(1) 


(2) 


(3) 


(4) 


(5) 


Arp 220E 


3.5 


2.2 




85 


Arp 220W 


7.5 


4.8 




95 


Arp 220E+W 


11.0 


6.9 


0.01 


95 


Mrk 231 


9.5 


33.5 


0.03 


1.5 


IRAS 08572+3915NW 


<4.0 


<30.0 


<0.05 


<3 


VV 114E 


3.0 


2.5 




20 


VV 114Esw 


2.5 


2.0 




110 



Note. — Col.(l): Object name. Col. (2): Observed flux of the 3.3 /im PAH emission. Col. (3): Observed 
luminosity of the 3.3 /im PAH emission. Col. (4): Observed 3.3 fim PAH-to-infrar ed luminosity ratio in units 
of 10 ~^ , a typical value for a modestly-obscured (Ay < 15 mag) normal starburst llMouri et al.|[l99ol : llmanishil 
l2002h . For Arp 220E, W and VV 114E, Eswi no values are shown, because the fraction of /Rj45'-measured 
infrared fluxes, coming from individual nuclei, is unknown. C ol. (5): Rest-frame equivale nt width of the 3.3 /im 
PAH emission. Normal starbursts have EW3.3PAH ^ 100 nm lllmanishi fc DudlevSoOOl) . 
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Fig. 1.— Continuum maps of Arp 220 at ~87 GHz and Mrk 231 at ~85 GHz. The 
crosses show the coordinates of main nuclei. The coordinates in J2000 are (15''34™57.22*, 
+23°30'11'.'5) for Arp 220W, (15'^34™57.29^ +23°30'11'.'3) for Arp 220E, and (12'^56™14.22^ 
+56°52'25'.'l) for Mrk 231. Thes e nuclear coordinates are adopted from C0(2-l) peaks in 
millimeter interferometric maps (IDownes fc Solomon! 119981 ). The contours are 3 x (3, 5, 7, 
9) mJy beam-i for Arp 220, 3 x (3, 5, 7, 9, 11) mJy beam"^ for Mrk 231. 
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Fig. 2 — Le/t: HCN(l-O) emission map. The contours are 0.85 x (3, 7, 11, 15, 19, 23, 27) Jy km 
s"^ beam~i for Arp 220, 0.76 x (4, 6, 8, 10) Jy km s"^ beam'^ for Mrk 231, 0.6 x (3, 4, 5) Jy km s"^ 
beam-i for IRAS 08572+3915, 0.5 x (3, 4, 5, 6, 7, 8) Jy km s'^ beam"! for VV 114, and 1.2 x (3, 
4, 5) Jy km s~^ beam~^ for He 2-10. For IRAS 08572+3915, the two crosse s correspond to th e NW 
(9''0"'25.38^ +39°3'53'.'8) and SE (9''0™25.62^ +39°3'48'.'8) peaks in J2000 (|Evans et aklbood l. For 
VV 114, the two crosses are near-infrared peaks, VV 114E (1^7™47.48", -17°30'25'.'6) andVVlMW 



(rT"46.53^ -17°30'22'.'l) in J2000, whose coordinates are converted from B1950 (IKnop et al 



1^994), using NED. The open star is CO(3-2) peak (l''7'"47.5^ -17°30'25'.'0) in J2000 (jlono et al 



2004 ). and the open tri angle is a radio 3.55 cm (8.44 GHz) peak (1^7™47.58^ -17°30'25'.'6) in J2000, 

using NED. E-1 and E-2 correspond to the two emission 



converted from B1950 (Condon et al 



1991 



peaks in the HCN(l-O) map. Fo r He 2-10, the cross is C O (1-0) peak (8^36^15. l^ -26°24'34'.'0) 



in J2000, converted from B1950 (jKobulnickv et al.lll995l ) using NED. Right: HCO+(1-0) emission 
map. The contours are 0.85 x (3, 5, 7, 9, 11) Jy km s"^ beam^^ for Arp 220, 0.76 x (3, 4, 5, 
6) Jy km s"^ beam'^ for Mrk 231, 0.6 x 3 Jy km s"^ beam"! for IRAS 08572+3915, 0.5 x (3, 
4, 5, 6, 7, 8, 9, 10, 11, 12) Jy km s'^ beam"! for VV 114, and 1.2 x (3, 4, 5, 6, 7, 8, 9) Jy km 
s~^ beam~^ for He 2-10. For VV 114, the positions 3 and 4 correspond to the two emission peaks 
in the HCO"''(1-0) map. For IRAS 08572+3915 and He 2-10, possible continuum emission is not 
subtracted from both the HCN(l-O) and HCO+(1-0) maps (see §4.1). 
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Fig. 3.— Channel maps of HCN(l-O) and HCO+(1-0) emission for Arp 220. Left: HCN(1- 
0) emission. The contours are 4 x (—3, 3, 5, 7, 9, 11) mJy beam~^. The r.m.s. noise level 
is ~4 mJy beam^^. Right HCO'^(l-O) emission. The contours are 4 x (—3, 3, 4, 5, 6) mJy 
beam"-*^. The r.m.s. noise level is ~4 mJy beam~^. 
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Fig. 4.— Channel maps of HCN(l-O) and HCO+(1-0) emission for VV 114. Left: HCN(1- 
0) emission. The contours are 3 x (—3, 3, 4, 5) mJy beam~^. The r.m.s. noise level is ~3 
mJy beam"^. Right: HCO+(1-0) emission. The contours are 3 x (—3, 3, 4, 5, 6, 7) mJy 
beam""^. The r.m.s. noise level is ~3 mJy beam~^. 
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Fig. 5. — HCN(l-O) and HCO"'"(1-0) spectra of the observed LIRGs. The abscissa is the 
observed frequency in GHz and the ordinate is flux in mJy beam^^. For Arp 220 and Mrk 
231, the strong continuum emission is subtracted. For VV 114, a possible continuum is 
subtracted. For IRAS 08572+3915, the continuum is not subtracted, because its subtraction 
introduces a significant uncertainty, given the faintness of the HCN(l-O) and HCO+(1-0) 
emission lines. Appropriate spectral binning is made, depending on the signal levels of the 
HCN(l-O) and HCO+(1-0) emission hnes. The r.m.s. noise levels per spectral bin in the 
final spectra are ~4 mJy (Arp 220), ~6 mJy (Mrk 231), ~1.5 mJy (IRAS 08572+3915), and 
~3 mJy (VV 114). 
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Fig. 6. — Gaussian fits to the detected HCN(l-O) and HCO+(1-0) emission lines. The 
abscissa is the LSR velocity {vopt = — 1) x c} in km s^^ and the ordinate is flux in 
mJy beam~^. Single Gaussian flts are used as defaults and are shown as dashed lines. The 
spectrum of IRAS 08572+3915 is only for the NW nucleus. For this source, since evidence for 
double peaks exists, we attempted a double Gaussian fit, shown as dashed fines. A constant 
continuum is assumed and is set as a free parameter because the continuum emission was 
not subtracted from the spectrum of IRAS 08572+3915. For other sources, the continuum 
level is set as the zero value. The adopted continuum levels are shown as the horizontal solid 
straight lines for all sources. 
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Fig. 7.— Newly obtained infrared 2.8-4.1 /im (L-band) spectra of Arp 220, Mrk 231, IRAS 
08572+3915NW, and W 114. The abscissa is observed wavelength in //m, and the ordinate 
is Fx in 10~^^ W /xm"-*^. Dashed straight line: Adopted continuum level to measure 
the optical depth of the 3.1 ^m ice absorption feature. Dotted line: Gaussian fit to the 
3.3 fim PAH emission feature. Solid straight line: Adopted continuum level to measure the 
strength of the 3.3 fim PAH emission feature. Dotted straight line for IRAS 08572+3915 is 
the adopted continuum level to measure the optical depth of the 3.4 /im carbonaceous dust 
absorption feature. 
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Fig. 8.— HCN(1-0)/HCO+(1-0) (ordinate) and HCN(1-0)/CO(1-0) (abscissa) ratios in bright- 
ness temperature (oc x flux density), derived from our NMA/RAINBOW interferometric ob- 
servations. Arp 220, Mrk 231, IRAS 08572-^3915, VV 114, and He 2-10 a r e plo t ted as large 
filled stars with labels. Oth er LIRGs previously observed by llmanishi et al.l (|2004l . l2006bl ) and 



Imanishi Nakanishil (j2006l ) are also plotted as small filled stars. Other data points are taken 
from lKohnJ ^20od ). where sources with AGN-like (starburst-like) ratios are marked with filled 
squares (open circles). The rat io at the core of the s tarburst galaxy NGC 253, derived by spatially- 
resolved interferometric data (IKnudsen et al.l 120071 ) . is also plotted as a large open triangle. For 
LIRGs and He 2-10, the HCN(1-0)/HCO''~(1-0) brightness-temperature ratios are derived directly 
from our interferometric data, while the HCN(1-0)/CO(1-0) ratios are derived by combining our 
data with CO data in the literature. The HCN(1-0)/CO(1-0) ratio in Arp 220 is taken as an 
upper limit, because the beam size of CO observations is much smaller. For all LIRGs' nuclei, the 
HCN(1-0)/HCO"'"(1-0) brightness-temperature ratios in the ordinate are those toward the cores, 
where putative buried AGNs are expected to reside in. Contamination from extended starburst 
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Fig. 9. — Schematic diagram of molecular gas in the current standard model (jSolomon et al. 



19871 ). Molecular gas consists of dense clumps (shown as dark filled circles in the upper-left 
panel) with a small volume filli ng factor. HCN(l-O) and H CQ"'"(l-0) lines can be moderately 
optically thick for each clump (iNguyen-Q-Rieu et al.lll992n . It is very likely that each clump 
has a decreasing ra dial density profile (oc r"; a < 0). If the radial density distribution 
of r~^'^ is assumed (iGierens et al.lll992l ). line em ission comes mainly from the layer where 
an optical depth reaches about unity (r ~ 1) (IGierens et al.l Il992l ). When a molecular 
abundance increases (decreases), the r ~ 1 layer for certain molecular line moves outward 
(inward). Consequently, the surface area of the layer, or its area filling factor in molecular 
gas, becomes larger (smaller), increasing (decreasing) the flux of that molecular line. 
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Fig. 10. — Comparison between FWHM line widths of molecular lines (abscissa) and 
HCN(1-0)/HCO+(1-0) brightness-temperature ratio (ordinate). Filled stars: HCN(l-O) line 
width. Open triangles: HCO+(1-0) line width. All LIRGs with available N MA/RAINBOW 



HCN(l-O) and HCO+(1-0) data are plotted. Th ese LIRGs are NGC 4 418 (llmanishi et al. 



2nn4h. UGC 5101, Mrk 273 IRAS 17208-0014 fllmanishi et al.ll2006bh . Arp 299 A, B, C 
(llmanishi fc Nakanishil 120061 ). and the four LIRGs (This paper). The HCN(l-O) line width 
of UGC 5101, and the HCN(l-O) and HCO+(1-0) line widths of IRAS 08572+3915 are for 
individual components of the double Gaussian fits (Imanishi et al. 2006b; This paper). The 
HCN(1-0)/IICO^(1-0) brightness-temperature ratios are also derived for individual compo- 
nents of the double Gaussian fits (Imanishi et al. 2006b; This paper). For Mrk 273 and 
VV 114E-1, no HCO^(l-O) line width information is available (Imanishi et al. 2006b; This 
paper). The FWHM line widths of NGC 4418 are 250 km s"^ [HCN(l-O )] and 170 km s"^ 
[HCO'^(l-O)], which are estimated from the data of llmanishi et al.l (120041 ). 




Fig. 11. — Comparison between /i?^5-measured infrared luminosity in units of the solar 
luminosity (abscissa) and HCN(1-0)/HCO+(1-0) brightness-temperature ratio (ordinate). 
The HCN(1-0)/HCO^(1-0) ratio in the ordinate is only for the core emission within the 
beam size of interferometric maps, while the infrared luminosity in the abscissa traces the 
emission from the whole galaxy regions. For VV 114, we plot the HCN(1-0)/HCO"'"(1- 
0) brightness-temperature ratio at VV 114E-1, and assume that the infrared luminosity of 
VV 114 is dominated by VV 114E. For UGC 5101 and IRAS 08572+3915, the HCN(1- 
0)/HCO"'"(1-0) brightness-temperature ratios are derived by summing the both components 
of the double Gaussian fits. For Arp 220, we assu me that Arp 220 W and E have, respectively, 
75% and 25% of the total infrared luminosity (IDownes fc EckartI I2007D . Fo r Arp 299, the 
infrared luminosities of Arp 299 A, B, C are taken from ICharmandaris et al.l (l2002l ). 



